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Abstract 

January 2009 saw the successful launch of the first space-based mission specifically designed 
for measuring greenhouse gases, the Japanese Greenhouse gases Observing SATellite (GOSAT). 
We present global land maps (Level 3 data) of column-averaged CO 2 concentrations (Xco 2 ) 
derived using observations from the GOSAT ACOS retrieval algorithm, for July through 
December 2009. The applied geostatistical mapping approach makes it possible to generate 
maps at high spatial and temporal resolutions that include uncertainty measures and that are 
derived directly from the Level 2 observations, without invoking an atmospheric transport model 
or estimates of CO 2 uptake and emissions. As such, they are particularly well suited for 
comparison studies. Results show that the Level 3 maps for July to December 2009 on a 1° x 
1 .25° grid, at six-day resolution capture much of the synoptic scale and regional variability of 
Xco 2 , in addition to its overall seasonality. The uncertainty estimates, which reflect local data 
coverage, Xco 2 variability, and retrieval errors, indicate that the Southern latitudes are relatively 
well-constrained, while the Sahara Desert and the high Northern latitudes are weakly- 
constrained. A probabilistic comparison to the PCTM/GEOS-5/CASA-GFED model reveals that 
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the most statistically significant discrepancies occur in South America in July and August, and 
central Asia in September to December. While still preliminary, these results illustrate the 
usefulness of a high spatiotemporal resolution, data-driven Level 3 data product for direct 
interpretation and comparison of satellite observations of highly dynamic parameters such as 
atmospheric CO 2 . 

1 Introduction 

The Greenhouse gases Observing SATellite “Ibuki” (GOSAT) launched on January 23, 

2009, and is the first space-based mission to reach orbit that was designed specifically for 
making high-precision measurements of carbon dioxide (CO 2 ) and methane (CH4) with 
sensitivity in the lower troposphere [Kuze et al. 2009; Yokota et al. 2009]. After the launch 
failure of the Orbiting Carbon Observatory (OCO) mission (e.g. Crisp et al. 2004), the OCO 
team was invited to join the GOSAT team in analyzing GOSAT observations, under the auspices 
of the NASA Atmospheric CO 2 Observations from Space (ACOS) task. The ACOS GOSAT 
column CO 2 (Xco 2 ) retrieval algorithm has now reached a level of maturity that makes it possible 
to use its estimates for informing carbon cycle science [O’Dell et al. 2012; Crisp et al. 2012], 
Version 2.9 of the Level 2 data product, which represents geo-referenced Xco 2 observations, 
includes approximately 900 successful retrievals per three-day repeat cycle during the second 
half of 2009, the first period for which data are available. The majority of these observations are 
over land. 

Although these data are useful in their own right, they have large gaps (e.g. Figure la and 
Figure Al in the auxiliary materials) and substantial retrieval uncertainties [O’Dell et al. 2012], 
which makes it difficult to interpret their scientific significance without further analysis. 
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Hammerling et al. [2012] recently developed a statistical mapping approach that makes it 
possible to create full-coverage (i.e. Level 3) maps from satellite Xco 2 observations at high 
spatial and temporal resolutions. Unlike commonly used spatial and temporal binning and 
averaging procedures (e.g., Crevoisier et al. [2009]; Kulawik et al. [2010]; Tiwari et al. [2006]), 
this approach exploits the spatial correlation among the Level 2 observations and the resulting 
Level 3 product describes the Xco 2 concentrations as a stochastic field characterized by its mean 
("Level 3 estimates") and variance ("Level 3 uncertainties") structure. 

Furthermore, unlike maps derived from inverse modeling or data assimilation studies (e.g. 
Engelen et al. [2009]), the Hammerling et al. [2012] approach draws information about the 
degree of spatial variability of Xco 2 directly from the Xco 2 observations, without additional 
information introduced from an atmospheric transport model or CO 2 flux estimates. As such, 
because no information from atmospheric transport models or CO 2 flux estimates is 
incorporated, the resulting Level 3 maps are a more direct representation of the information 
content of the retrievals. Rather than being intended as inputs to inverse modeling studies, these 
Level 3 Xco 2 products enable direct independent comparisons with existing models of carbon 
flux and atmospheric transport. The uncertainty measures provided by the approach make it 
possible to conduct these comparisons in a probabilistic framework. 

This paper presents global Level 3 X C 02 products over land derived from the GOSAT ACOS 
Xco 2 retrievals, covering the second half of 2009. The Level 3 estimates and their associated 
uncertainties are compared to predictions for the same period from a combined CO 2 flux and 
atmospheric transport model using a probabilistic framework. 

2 Data and Methods 
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2.1 GOSAT ACOS X C0 2 Level 2 data 


GOSAT flies in a sun-synchronous orbit with an approximate 1pm equator-crossing time and has 
a three-day repeat-cycle. Version 2.9 of the GOSAT ACOS Xco 2 Level 2 data product is used in 
this study; only high (H) gain data were used as recommended in Crisp et al. [2012], Figure la 
shows an example of six days (i.e. two repeat cycles) of ACOS L2 data for August 2009. 

2.2 Method for creating global GOSAT ACOS Xco 2 Level 3 maps 

The geostatistical methodology applied for creating Level 3 maps exploits the spatial 
correlation of the Xco 2 observations and consists of two major steps. In the first step, the spatial 
covariance structure of the Xco 2 observations is inferred from these observations. In the second 
step, the inferred spatial covariance structure and the observations are used to estimate the Xco 2 
field. The approach is described in detail by Hammerling et al. [2012], and only key 
implementation details are presented here. Due to the currently limited availability of GOSAT 
ACOS Xco 2 observations over the oceans, the estimation has been restricted to land areas. The 
mapping is implemented on a 1° latitude x 1.25° longitude grid, to inform regional variability 
and to correspond with that of the model used for comparison in Section 4. 

Based on previous work [Alkhaled et al., 2008], an exponential covariance function is used 
to represent the Xco 2 spatial correlation: 

- ( h \ (1) 
C(f£) = o- exp | - jj, 

\ - ' 


where the covariance C is a function of the separation distance between locations (A), and 
spatially- variable variance {a 2 ) and range (/) parameters that are inferred at each estimation 

location from the Level 2 data. 
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A local kriging procedure is then applied to create full-coverage maps, using a weighted 
average of available observations by solving the following linear system of equations once for 
each location on the Level 3 map: 


Q-R 

l r 



( 2 ) 


where Q is an n X n covariance matrix among the n observation locations, as defined in eq. 1, R 
is an n X n diagonal matrix with the retrieval error variance specific to each observation on the 


diagonal, A is a « X 1 vector of weights, is a Lagrange multiplier and q is the n X 1 vector of 

the spatial covariances between an individual estimation location and the observation locations, 
also defined using eq. 1. In this study, the measurement error variances are the squares of the 
reported ACOS Level 2 measurement error standard deviations adjusted by a factor of 2. 1 as 
derived by O'Dell et al. [2012]. The predicted Xco 2 value, y, and the prediction uncertainty, er 2 ^, 

at each Level 3 location are: 


A r y 

(3) 

= a 2 — Aq — v. 

(4) 


where y are the observations at the n Level 2 locations and a 2 is the variance as shown in eq 1 . 

Based on previous work, a 2000 km neighborhood is required for assessing the local 
spatial variability (eq. 1, also see Hammerling et al. [2012] for details), and estimates can 
therefore only be obtained if there is a minimum of three observations within this distance of 
each estimation location. Estimation locations not meeting this requirement are shown as white 
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in Figures lb and A2. It is the uncertainties in eq. 4, however, that should be used as the 
criterion for limiting the coverage of Level 3 maps to regions where they are interpretable for a 
given scientific application (e.g. Figure 2), and one of the advantages of the method is the 
flexibility to dynamically define this uncertainty tolerance. 

2.3 PCTM/GEOS-5/CASA-GFED model data 

The modeled Xco 2 data used in the intercomparison are based on the Goddard Space Flight 
Center parameterized chemistry and transport model, which is driven by real-time analyzed 
meteorological fields from the Goddard Global Modeling and Assimilation Office, version 
GEOS-5, and uses biospheric fluxes produced from the Camegie-Ames-Stanford-Approach, 
which incorporate biomass burning from the Global Fire Emissions Database (PCTM/GEOS- 
5/CASA-GFED), as well as oceanic, and anthropogenic CO 2 flux estimates, as described by 
Kawa et al. [2004; 2010], The model resolution is 1° x 1.25° with 28 vertical levels and hourly 
output. CO 2 mixing ratios were pressure-averaged to simulate the vertical sensitivity of the 
GOSAT observations. The PCTM/GEOS-5/CASA-GFED model has been widely tested, and 
has shown favorable results in carbon cycle comparison studies (e.g. Kawa et al. [2010] and 
references therein). 

3 ACOS GOSAT Xco 2 Level 3 maps 

The choice of the temporal resolution, meaning the time period over which observations are 
aggregated, is an important decision in the creation of a Level 3 product [Hammerling et al., 
2012], Ideally Level 3 products are created for the shortest time period possible to preserve as 
much of the short-term dynamical information as possible. However, this needs to be balanced 
with a minimum requirement for spatial coverage by the GOSAT observations. Based on initial 
investigations of temporal resolutions ranging from three days to one month, a resolution of six- 
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days fulfilled both these objectives for all 30 six-day periods investigated from July to December 
2009. 

Figure 1 provides an example of one of the investigated periods, August 7-12 2009. The 
Level 3 map (Figure lb) for this period shows comparatively low Xco 2 in the Northern latitudes 
consistent with the knowledge of the effect of the seasonal cycle on CO 2 concentrations. The 
fact that the seasonal cycle in the Northern hemisphere is captured well in the GOSAT ACOS 
Level 3 maps becomes further evident from results from the frill examined period (see Figure A2 
in auxiliary materials), which show a pronounced increase in CO 2 concentrations in the Northern 
latitudes in the winter months as well as a more subtle increase in the overall C0 2 
concentrations. The comparatively high Xco 2 over South America visible in the Level 3 map for 
August 7-12 (Figure lb) are a fairly persistent feature throughout the summer months (figure A2) 
and are further discussed in Section 4. 

An advantage of the mapping method used in this study is that each estimate has an 
associated uncertainty measure (Figure lc), which reflects the number of observations 
surrounding an estimation location, their retrieval errors, and the spatial variability in the Xco 2 
field. Locations where the prediction uncertainties are below specific cut-off values are 
illustrated in Figure 2. For this six-day period, the predictions uncertainties are low for 
Australia, the southern part of Africa and eastern South America, whereas they are high for 
Southeast Asia, parts of India and the eastern United States and Canada. Analyzing these 
prediction uncertainties over extended time periods highlights the degree to which ACOS 
GOSAT retrievals constrain the Xco 2 distribution for different regions. Figure 3 summarizes this 
analysis for the 30 investigated six-day periods in 2009, identifying Australia, Southern Africa 
and a region in South America covering approximately eastern Brazil, Paraguay, Uruguay and 
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central and northern Argentina and northern Chile as well-observed regions. Regions with the 
weakest constraint are the Sahara Desert and the high Northern Latitudes including Alaska, 
northern Canada, Greenland, Scandinavia and northern Russia. 

The interplay of how the number of observations, their retrieval errors, and the spatial 
variability in the Xco 2 field contribute to the uncertainty at each location renders it difficult to 
completely separate the effect of these contributing factors. The spatial coverage over the larger 
land masses in the Southern hemisphere, namely Australia, southern Africa and southern South 
America, is generally good. The number of observations decreases somewhat towards the end of 
the year, but these observations have lower retrieval errors and are supplemented by nearby 
ocean observations, which shift southwards in the second half of year as a function of the solar 
zenith angle. Southeast Asia and central and eastern China, on the other hand, have very poor 
coverage during July to October due to persistent cloudiness, but notably better coverage in 
November and December, leading to the mapping uncertainties being seasonally variable. The 
United States have generally good coverage, but the Xco 2 spatial variability over the Northern 
hemisphere land masses is rather high, yielding somewhat higher uncertainties for North 
America than for areas with comparable spatial coverage but less spatial variability such as 
Australia. There are no observations over the Sahara Desert, due to our exclusion of the GOSAT 
M-gain data (see Crisp et al. [2012] for details). The high Northern latitudes lack observations in 
November and December due to solar zenith angle restrictions; and the observations in July to 
October have comparatively high retrieval errors. This, coupled with the high Xco 2 spatial 
variability in the high Northern latitudes, leads to high mapping uncertainties even when data are 
present. 
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4 Comparison of Level 3 maps to modeled X C o 2 

The ACOS GOSAT Xco 2 Level 3 products can be used to conduct intercomparisons with 
models, by using the Level 3 data and their associated uncertainties to probabilistically identify 
areas where model outputs differ significantly from the Level 3 maps. 

Figure 4 shows an example of such an intercomparison to the PCTM/GEOS-5/CASA-GFED 
model for August 7-12 2009. The difference plot (Figure 4b) shows large differences in North 
America, the Amazon Region, and in a region covering the Northeastern part of India and 
Bangladesh. The standardized differences (Figure 4c), on the other hand, incorporate the Level 3 
uncertainties, and can therefore be used to assess the significance of these differences given the 
information content of the satellite observations. For example, while the difference in North 
America and Southeast Asia might appear large in Figure 4b, they are not highly significant, as 
shown in Figure 4c. This is due to the comparatively large Level 3 uncertainties in these regions 
for this period. 

Figure 5 summarizes the intercomparison for July to December 2009, and reveals that 
discrepancies are most pronounced over South America for the Northern hemisphere summer 
months and shift to Asia in the Northern hemisphere fall. Although these results likely point to 
areas where the PCTM/GEOS-5/CASA-GFED model flux and transport processes need to be re- 
examined, Level 2 retrieval biases and, in the case of the sparsely-sampled Amazon region, 
underestimation of the Level 3 uncertainties due to low Xco 2 variability in surrounding well- 
sampled regions cannot be absolutely eliminated at this stage. It is also interesting to note that 
certain regions exhibit few or no limited significant differences over the entire examined period, 
including the high Northern latitudes, North America, Northern Africa, the Arabian Peninsula 
and Australia. The conclusion one can draw from an absence of statistically significant 
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discrepancies depends on how well constrained a region is. For example, the high Northern 
latitudes are weakly constrained and have high mapping uncertainties. This implies that even 
large discrepancies are not conclusive because the power to detect a difference is low for that 
region. For Australia, on the other hand, the Level 3 uncertainties are rather low, so an absence 
of detectable discrepancies indicates that the Level 3 maps are indeed consistent with the model 
outputs. 

5 Conclusions 

This paper presents global Xco 2 Level 3 products over land based on the ACOS GOSAT 
Xco 2 data. The implemented approach (Hammerling et al., 2012) yields maps at high spatial and 
temporal resolutions, using information derived directly from the Level 2 observations, without 
invoking an atmospheric transport model or estimates of CO 2 uptake and emissions. One 
limitation of such a purely observation-driven approach is that local enhancement phenomena 
that are not observed by the satellite cannot be fully captured. This results in Level 3 maps with 
smoother features than expected in the real Xco 2 concentration fields, but with uncertainty 
bounds that are wide enough to capture the range of likely variability. 

Level 3 maps for July to December 2009 at six-day resolution capture much of the synoptic 
scale and regional variability of Xco 2 , in addition to the overall seasonality. Results include 
robust uncertainty estimates, which reflect local data coverage, Xco 2 variability, and retrieval 
errors. Uncertainties are generally highest in the northern hemisphere in July and August, during 
the height of the growing season (Figure A3), and lowest in areas with good data coverage and 
low CO 2 variability in the Southern Hemisphere (Figure A3). 

A probabilistic comparison to a state of the art model reveals that the most significant 
discrepancies captured by the ACOS GOSAT Level 3 maps are in the South America in July and 
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August, and central Asia in September to December (Figure 5). The differences in South 
America are significant in part because the Level 3 uncertainties are low in this region, and may 
reveal inaccuracies in carbon flux estimates for this region that is poorly constrained by in situ 
atmospheric CO 2 observations, although problems with the Level 2 retrievals and with 
identifying local phenomena in the Amazon region in the Level 3 products cannot be ruled out at 
this stage. Similarly, the significant differences in Asia appear during months when the Level 3 
mapping uncertainties are lowest in this region (Figure A3). 

These early results illustrate the usefulness of a high spatiotemporal resolution, data-driven 
Level 3 data product with uncertainty measures. Such a Level 3 data product can be used for 
direct interpretation of satellite observations, including those of highly dynamic parameters such 
as atmospheric CO 2 , and for probabilistic comparison studies. 
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Figure Captions 

Figure 1: ACOS Xco 2 Level 2 data ("Observations"), ACOS Xco 2 Level 3 product 
(“Estimates”) and estimated prediction uncertainties (“Uncertainty”) expressed as a standard 
deviation for August 7-12 2009. 

Figure 2: ACOS Level 3 X C o 2 map for August 7-12 (Figure lb) filtered for locations where 
the standard deviations of the prediction uncertainties (Figure lc) are less than 2.5ppm, less than 
2ppm and less than lppm, respectively. 

Figure 3: Summary of the analysis of prediction uncertainties from 30 6-day periods from 
July through December 2009. For each location, the number of 6-day prediction periods with 
prediction uncertainties below 2.5 ppm, 1.5 ppm and 1 ppm, respectively, is shown. Lighter 
colors indicated regions which are better constrained by the GOSAT observations. 

Figure 4: PCTM model predictions for the same 6-day period as shown in Figure 1, 
difference and discretized standardized difference between the ACOS Level 3 map and the 
PCTM model. In the difference plot, values in the copper range indicate areas where the ACOS 
Level 3 values exceed the PCTM model predictions, values in the blue range areas where the 
PCTM model exceed the ACOS Level 3 values. The standardized difference is the absolute 
difference divided by the standard deviation of the prediction uncertainty at each location. The 
values are discretized to improve the visualization. Areas in yellow represent differences larger 
than one standard deviation of the prediction uncertainty, areas in orange larger than two 
standard deviations and areas in dark red larger than three standard deviations. The PCTM data 
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has been mean-adjusted to the ACOS Level 3 predictions so that the global spatial average of the 
PCTM data and the ACOS Level 3 is equal. 

Figure 5: Percentage of 6-day periods within each month where the standardized differences 
exceed two prediction uncertainties. Given that there are only five 6-day periods in each month, 
we have chosen a discretized color scale to visualize percentages. 
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